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Abstract 


A barotroplc model of the atmosphere Is used to test various sources 
of forecast error. These errors are classified as truncation error, 
physical error, or Initial error. It Is shown that growth patterns 
due to each category differ significantly. Initial er-ors are sho’in 
not to grow In a barotroplc model contrary to reports of other studies 
which Indicate that they basically do grow. Also, random Initial 
errors are shown to decrease due to the filtering effect of the model 
Itself, Results seem to Indicate that Instabilities are required for 
error growth, be they barotroplc or barocllnlc, and that random errors 
are not representative of true Initial conditions. 


INTRODUCTION 


With the launching of Soasat In Spring, 1978, a significant amount 
of meteorological data will bo added to tho forecasters Inventory. Two 
of tho satellite's sensors will provide much needed Information on tho 
oceans' surface temperature and wind stress. These may bo Incorporated 
Into tho Initial data sets of large numerical forecasting models and 
hopefully reduce errors In forecasting. Before estimates con be made, 
however, on the effectiveness of these new data, or of any satellite 
observations, for that matter, It Is necessary to understand tho nature 
of forecasting errors, In general, and the role Increased observations 
can ploy In reducing them. 

Errors In weather forecasts are generally broken down Into three 
main categories: 

1. Physical error due to poor representation of the actual 
phenomena by the mathematical equations. 

2. Truncation error due to finite difference approximation of 
the governing equations, and 

3. Initialization error due to Inaccurate portrayal of the 
Initial conditions. 

The physical error Is generally due to a lock of knowledge concern- 
ing the physical processes of the atmosphere, or, because of computer 
limitations, certain phenomena aro compromized with coarse porametorlza- 
tlons. That Is to say, that processes, such as radiation, for example, 
may require frequent calculation of small components such as absorption 
by wavelength bands, but, because of time limitations, the wavelengths 
are broken Into only 3 or 4 bands and tho absorption and transmission at 
these bands treated with some simple, bulk coefficients. Because these 
parameter {nation s are, by definition, only approximations to the 


true case, one would expect some error to result from their use. How 
great this error may bocome con only be approximated by simulations 
where more precise equations are used for a limited time. In processes 
where no exact mathematical representation Is known, no real assessment 
of the physical error can be furnished. 

Truncation errors are a general category embodying all errors that 
may result from the use of a coarse grid for a finite difference solution 
to the governing equations. In a mathematical sense, the errors are due 
solely to the use of finite-differences Instead of exact derivatives. 

From a physical perspective, the use of a coarse grid filters many of 
the small scale phenomena which may still be Important to weather forecasting. 
Mlyakoda, et al., (1971) have demonstrated the effects of Increasing model 
resolution by forecasting with several versions of the Geophysical Fluid 
Dynamics Laboratory (GFDL) 9-level general circulation model (GCM) . 

There are significant Impacts on the forecast when the grid spacing Is 
decreased from 5^0 to 270 km but much smaller Impacts when ft Is further 
decreased to 135 km. Apparently the processes which are omitted In the 
very coarse grid are quite Important even for short-range forecasts, 
while there Is little contribution to the short-range forecasts from the 
smaller scale. The 135 km-grld did become significant after about 5 to 
6 days, however, proving that even small scale phenomena are Important 
for longer range forecasts. 

Initial errors Include all deviations of the model's Initial state 
from the true Initial state. These deviations can be caused by random 
Instrument error, human error, omission of large areas, and the sub- 
sequent Interpolation errors Involved In applying the sparse observa- 
tions to the grid points. These various sources may not all result In 
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the same typo of error distribution. In fact, the random Instrument 
error Is generally of low magnitude and Independent of location. The 
other errors (excluding human error) ore geographically dependent, 
with the oceans and uninhabited areas notoriously lacking In meteorological 
observations. 

In practice, It Is difficult to separate the Influences of these 
Individual eri\>r sources, although attempts have been made to simulate 
the various error patterns associated with each. Much of the concern 
over error growth has been stimulated by the necessity to Incorporate 
satellite data Into numerical models of the atmosphere; Kasahara (1972) 
presents a summary of various studies on the nature of data assimila- 
tion Into numerical models and briefly reviews the progress In under- 
standing Initial error growth. Bongtsson (1975), too, discusses the 
nature of errors In relation to the assimilation of satellite data Into 
numerical models of the atmosphere as does Blumen (1976a) and Blumen 
(1976b). The work here Intends to compar>. the various sources of 
error and their effects on a simplified, non-! Inear model of the atmos- 
phere. Understanding these various errors will help In the assessment 
of Impacts on forecasts due to the assimilation of current and projected 
satellite data. Part of this study was assimilated Into the Seasat simula- 
tion experiments conducted at the Goddard Institute for Space Studies 
(GISS) and described by Cane, et al. (1977). 

II. DESCRIPTION OF THE EXPERIMENT 

The model selected for studies of error growth was the simple, 
shallow water equations, which assumes a barotroplc, Invlscld, Incom- 
pressible fluid with height-independent velocities. These equations 
have been used for assimilation studies, In either their primitive 
form or In their differentiated vortlclty form, by many researchers 
Including Blumen (1976a), Morel, et a!., (1971), and Mlyakoda and Talagrand 
(1971), and Halberstam (197*1). The governing equations are: 


0) a, 3u/3t + u 3u/3x + v 3u/3y = -g 3h/3x - fv 

b, 3v/3t + u 9v/3x + v 3v/3y ® -g 3b/3y + fu 

c. 9h/3t + u 3h/3x + v 3h/3y » -h (3u/3x + 3v/3y) , 

where u and v are the westerly and southerly wind components, h Is the 
height of the fluid, f Is a corlolls parameter, and g Is earth’s (con- 
stant) gravitational acceleration. The boundaries were cyclical In the 
x-dlrectlon and Impenetrable, yet frlctlonle&s at y»0 and y°Y. f was 
assumed equal to 20 « 7.2 x 10 ^ s. The latitudinal variation In 
the model was from 20N to 70N. 

Most runs were made with a resolution of about 2-1/2° In the 
meridional direction and about 5° of longitude In the latitudinal direc- 
tion. This employed a grid system of 42 x 21 points. The numerical 
scheme featured a staggered grid (the "C n scheme of Arakawa and Mlntz, 
1974) with u, v and h on separate grids. The North and South boundaries 
were v-polnts with v set to 0, permanently. A centered difference time 
step equal to 10 min. was used, except once every 2 hours, when a Matsuno 
predictor-corrector step was used, to Insure against separation Into two 
computational modes. 

The simulated "perfect" run was made from prescribed Initial 
fields of h, u, and v, where h varied sinusoidally, while u and v 
were balanced geostrophlcally. The equations were Integrated for 240 
hours (10 days) and the fields stored every 2 hours. This run (hereto- 
fore referred to as "nature") was used as the verification for all 
other runs. The Initial distribution of h and the resultant 6-hour 
field are shown In Figure 1. 
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The "error" runs were* designed to test the relative magnitude and 


orror growth patterns associated with each type of error. The first 
of these was the "Initial" error Including all effects of wrongly pre- 
scribing the Initial field. Many of the simulation experiments men- 
tioned previously chose to perturb the Initial field by the addition 
of random errors. But many of the errors prevalent In the objective 
analysis of the numerical prediction models are not necessarily random, 
as mentioned earlier. The lack of data over wide expanses of the globe 
and the Interpolation and assimilation of rare bits of Information over 
large domains can easily be a more serious source of error than random 
Instrument or human error, for this reason, both types of Initial 
error wore tested separately here. The first type was generated by 
adding random errors to h, u, and v and allowing the modol to predict 
forward without further disturbance. This run will be designated as 
RE. The second {called DG for "data gap"), was generated by omitting 
Information from the center of the field and extrapolating linearly 
from the surrounding points, figs. 2 and 3 show the height con- 
tours at initial time and at 6 hours for RE and 06, respectively. 

Those should be compared to fig. 1 to assess the degree of departure 
from nature. Note that the Initial field for RE Is badly mottled, 
but some recovery Is noted by 6 hours. DG has an apparent "hole" In 
the middle of the Initial field which seems to cause a distortion Vn 
the field at 6 hours. 

The second type of orror Investigated was physical error, which 
can bo simulated In numerous ways. In the real world, It ts difficult, 
If not Impossible to evaluate the true physical error associated with 
Incorrect physical assumptions or mathematical parameterization. In 
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cases whore these con be tested, results ore usually too select to bo 
meaningful . Williamson and Kasahara (1970 hove attempted to show the 
offoet of physical error In relation to Initial error by demonstrating 
the effect of changing their boundary layer parameterl ration. But not all 
types of physical error behave In the same way. Obviously a ‘'dry" model 
does not produce the same forecast os one which Includes the effects of 
water vapor. In this experiment, wo tested two types of physical error. 
One Involved linearising the right-hand side of (l)c by substituting an 
average value of H for h. The other resulted from the Inclusion of a 
heating torm added to (l)c which effectively raised h In the lower part 
of the field and decreased It In the upper part. The Initial fields were 
the same as . >.'ire In both cases. The six hour h- field Is shown In Fig. 

I» for the linearised case (dubbod PE) and Is not much different from 
naturo; contours for the heated case ore not available. 

The third type of error Is a simple truncation error produced by 
decreasing the grid resolution from a *»2 x 21 network to \ 7 and 
Increasing the grid spacing threefold. The new "error" will thus 
represent a mathematical truncation relative to the fine resolution. 

No physical changes wore made In the model to compensate for the 
change In resolution and oven the time stop was not Increased, although 
linear stability considerations would have oMowed for It. The h 
contours for 0 and 6 hours are shown In Fig. 5 for the low resolution 
(IR) run. Note that Initially there Is considerable smoothing of the 
waves, with some of the closed contours appearing In nature missing In 
IR. At 6 hours, the gradients re >mble more the Initial conditions 
than the 6-hour naturo field, with the tight gradients still occupying the 
northern portion of the domain rather than moving down to the southern 
portion as found In nature. 



IN. ERROR GROWTH PATTERNS 

All tho error runs contl nuod until 2^0 hours and tho root moan 
square (rms) error for tho entire field was computed every 2 hours 
until 12 hours and every G hours thereafter, A plot of h-rms vs. time 
Is presented In Fig. G and u-rms In Fig.. 7 for PE, the heated case, 

LR, RE, and DG. They are obviously all quite dtfforent. 

Tho physical errors show a considerable range of possibilities. 

Tho PE run asymptotes at fairly moderate levels. After a rapid rise 
during tho first 6 hours, tho error growth rate declines until It bogtns 
to oscillate at about 2h to 36 hours. Apparently tho differences 
between tho fields tend to remain small but same waves may bo slightly 
out of phase. The heated cose, on the other hand, shows a continued 
growth of error for the entire 10 Joys, probably because the magnitudes 
of tho fields are considerably dtfforent In tho heated case. 

The LR case shows an almost predictable orror growth pattern. 

There Is an Inttla! spurt of error growth followed by a choppy os- 
cillation about some mean quantity. The oscillations are probably 
due to phasing differences caused by the differences In resolution 
between nature and LR. In linear wave equations, the truncation errors 
are also expected to oscillate since they are proportional to some order 
derivative of the time solution which consists of sinusoidal waves. 

The Initial errors present tho most Intriguing growth patterns In 
this study, RE shows a drop in error levels from the Initial time until 
some asymptote Is reached. The u-rms first Increases as It adjusts 
to the error In h, but drops along with h Immediately afterwards. DG, 
however, seems to maintain the same orror level with a possible In- 
crease towards the end of the period. 


One can understand the behavior of the Initial errors from a dis- 
cussion by Smagorinsky, et al., (1970)- who blame the growth of Initial 
error on barocltnlc Instabilities that amplify the Initial perturba- 
tions until the perturbed field resembles the original field only as would a 
random cl imataloglcal state . At the outset, the field goes through an 
adjustment period, according to Smagorinsky, et al., (1970), when the 
errors actually dip before rising. This dip was also mentioned by 
HaJberstam (197^) and Blumen (1976a) using the barotroplc equations. 

This study seems to Indicate that the "adjustment" Is a 
phenomenon associated only with random initial error. Other modes 
of error, similar to DG described here, may lead to different growth 
patterns where the noted dip does not occur. Apparently the dip 
is associated with the filtering of high-frequency waves by the numerical 
model. The noisy waves are a non-meteorologlcal component which con- 
taminate the field, and when they are removed by the smoothing and 
physics of the model, the field is brought into better agreement with 
other model fields or, at times, even with nature. 

The decrease of error with time in RE seems to support Smagor insky's, 
et al., (1970) contention that the barocllnlc instabilities are responsible 
for error growth. In our study, a barotropic model was used, so that the 
Initial errors either decrease or remain fairly constant. One must 
explain, however, why Miyakoda and Talagrand (1971) or Blumen (1976a) do 
show error growth even with barotropic models. They all instst that 
the non-linearity of the model is necessary for error growth, but It 
should not be a sufficient condition. In Miyakoda and Talagrand' s 
(1971) study, a barotropic instability may have been created by proper 


boundary conditions producing sufficient wind shear. Bl union's model 
separates the goostrophlc from the ageostrophlc components of tho 
vortlclty and his analysis shows that tho expected value of tho rms 
error, assuming a random Initial distribution, reachos an asymptotic 
valuo after a number of days, similar to realistic tests done with tho 
National Contor of Atmospheric Research (NCAR) model. Ills approach Is 
quite different from a stratgh forward numerical approximation to tho 
primitive equations or tho vortlclty form. His equations ore, In 
effect, linear, with non-llnoar components addod to them. This may have 
had an effect on the rate of computed error growth, although parallels 
to the primitive equations calculations should still exist. 

Morel, at al., (1970 also used the barotroplc case to study the 
offect of data assimilation on tho rate of error growth. Tho Initial 
conditions aro also porttrbod by random orrors and simulated satellite 
data assimilated during an Iterated forocastlng-hlndcastlng cycle. Tho 
error Is seen to decline during tho assimilation cyclo, but tho control 
caso (without the assimilation of data) Is omitted. Tho reader Is led 
to believe that tho docroasa In error Is strictly duo to tho assimila- 
tion of data, rather than to the adjustment process. In light of tho 
present study, their work may have to be reviewed with a non-asslml la- 
tod run presented In comparison. If the assimilations aro not totally 
responsible for the reduction of errors, then tho recommended, expen- 
sive Iterative process may not be worthwhile, especially with largo 
numerical atmospheric models. 

IV. SUMMARY AND CONCLUSIONS 

A primitive equations barotroplc modol of the atmosphere has boon 
used to tost the effects of certain errors on numerical forecasts. 

The throe typos of errors Investigated wore physical errors, truncation 
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errors* and Initial errors. The physical errors consisted of a change 
In the model equations by 1, substituting a constant average value of H 
for the variable h In the divergent term of 1 (c) and 2, Including a 
heating term which causes h to rise In the southern portion and de- 
crease In the northern portion. The truncation error was generated 
by reducing the resolution to 1/3 of the original l i2 x 21 grids. The 
Initial errors were divided Into a I. random error and 2. sampling and 
interpolation error created by omitting the center of the field and 
filling the gap by Interpolation. 

Results show that physical errors have a large range, depending 
on the type of physics Involved. Resolution errors grow Initially and 
then oscillate about some mean quantity. Initial errors do not grow 
In this barotroplc model. Instead, the random errors decrease until 
some asymptote Is reached, while the Interpolation error maintains 
a small oscillation about some mear value determined at the outset. 

Indications are that physical errors are difficult to assess and 
probably more difficult to correct. The ones that are due to Incorrect 
representation of the phenomena are the product of either unknown physics 
or unavoidable parameter I zat for! due to computer considerations. 

Truncation errors will always exist as long as ftnlte differences are 
employed. Mtyakoda, et a)., (1971) have clearly shown that there Is a 
critical resolution which must be surpassed In order fur the forecast 
to be meteorologically meaningful. Beyond that, very little incremen- 
tal benefit Is achieved by Increasing the resolution, at least for the 
short-term. 

Host efforts In Improving observations, such as the launching of 
satellites, are geared toward reducing the Initial errors. This study 
observes that not all initial errors behave In the same way, and that 
they will not grow In the absence of instabilities. Simulation studies 
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that make uso of only random Initial arrors, may not be doing justice to 
the problem of Initialization, slnco random errors are generally filtered 
out by the models themselves. On the other hand, areas and seasons of 
low barocl Inlcl ty may not require concentrated Improvement of Initial 
ftolds to Improve forecasts. 

Further Investigations ore planned to study the nature of error 
growth, especial ly wl th regards to Initial error. Experiments per- 
formed at the Goddard Institute for Space Studies (GlSS) with a 
realistic Initial Field are currently being documented by Cane, et al., 
(1977). A two-level harocltnlc model Is being developed at JPL to test 
the effects of barocl Inlcl ty on error growth, and whether areas of 
potentially serious errors can be predicted. 
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Fig. 7. Rms of u vs time between the nature field and (a) PE, 
(b) the heated case, (c) LR, (d) RE, and (e) DG 
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